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, Abstract. We consider the Cauchy problem for the Gerdjikov- 

Qh! Ivanov(GI) type of the derivative nonhnear Schrodinger (DNLS) 

"^ I equation: 

with stephke initial data: 17(2;, 0) = for a; < and q{x, 0) = 
j^^-2tBx Jqj. ^ ^ OjWhere A > and i? g M are constants. The 
paper aims at studying the long-time asymptotics of the solution 
to this problem. We show that there are four regions in the half- 
^ , 1 plane — cxD<a;<cx),i> 0, where the asymptotics has qualitatively 

different fornis:a slowly decaying self-similar wave of Zakharov- 
Manakov type for x > —AtB, a plane wave region:a; < —4t{B + 
, ^2A2(B+^)), an elliptic Tegion:-U{B+^2A'^{B + ^)) <x < 

\^ I —4tB. The main tool is the asymptotic analysis of an associated 



matrix Riemann-Hilbert problem. 



1. Introduction 

The classical, matheniatical model for non-linear pulse propagation 
in the picosecond time scale in the anomalous dispersion regime in an 
isotropic, homogeneous, lossless, non-amplifying, polarization-preserving 
single-mode optical fibre is the non-linear Schrodinger(NLS) equation 
[2]. However, in the subpicosecond-femtosecond time scale, experi- 
ments and theories on the propagation of high-power ultrashort pulses 
in long monomode optical fibres have shown that the NLS equation is 
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no longer valid and that additional non-linear terms (dispersive and 
dissipative) and higher-order linear dispersion should be taken into 
account, you can see [3S] and the references therein. In this case, 
subpicosecond-femtosecond pulse propagation is described (in dimen- 
sionless and normalized form) by the following non-linear evolution 
equation (NLEE) 

1 _ _ 7- 

iUf H Urr + \u\'^U + is{\u\'^u)r = —iVu + iSU-rTT H --u{\u\'^)ti (1-1) 

2 To 

where u is the slowly varying amplitude of the complex field enve- 
lope, ^ is the propagation distance along the fibre length, r is the time 
measured in a frame of reference moving with the pulse at the group 
velocity (the retarded frame), s(> 0) governs the effects due to the 
intensity dependence of the group velocity (self-steepening), T is the 
intrinsic fibre loss, 6 governs the effects of the third-order linear dis- 
persion, and — , where Tq is the normalized input pulsewidth and r„ is 
related to the slope of the Raman gain curve (assumed to vary linearly 
in the vicinity of the mean carrier frequency, Uq), governs the soliton 
self-frequency shift (SSFS) effect, [36] and the references therein. 

We set the right-hand side of (II. ip equal to zero, we obtain the 
following equation, 

m^ H — Urr + \u\'^u + isdwp-u),- = 0, (1.2) 

This equation is related to the Kaup- Newell type of derivative nonlinear 
Schrodinger equation, 

iqt{x, t) = -qxx{x, t) + {qq^)x (1-3) 

by change of variables 



u{,i, r) = q{x, t)e'^^ ^\ ^ = -^, r = -^ + -^. 



t X t 

2s"2' ^ ""2^^273 • 



And we note that if we replace x by —x, equation (II. 3p changes into 

iqt{x,t) = -q^^{x,t) - (qq^)^. (1.4) 
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But, we also know if we formulate a Rieinann-Hilbert problem for the 
solution of the inverse spectral problem of the equation fll.4p . we find 
we cannot find solutions of its spectral problem which approach the 
2x2 identity matrix I as A; — )■ cxD.It is well-known that there are three 
kinds of celebrated DNLS equations, including Kaup-Newell equation 
( i.e Eq.f ll.4p ). Chen-Lee-Liu equation [37j 

iqt + Qxx + ilql^Qx = 0, 

and Gerdjikov-Ivanov(GI) equation [381 SO] 

iqt + qxx-iq^qx + ^\q\'^q = (1.5) 

It has been found that they may be transformed into each other by 
gauge transformations [HHl |39]. And in |40], the Gl-type has the re- 
quired property of the solutions of its spectral problem which approach 
the 2x2 identity matrix I as A; — t- oo. So,we focus on the Gl-type of de- 
rivative nonlinear Schrodinger equation. In the following of this paper 
we also name the Gl-type DNLS equation as DNLS equation. 

Initial value problems for nonlinear evolution equations with step- 
like initial data have attracted much attention since the early 1970s 
[m [T71 [IHl [19], but only a few rigorous results concerning the long- 
time behavior of solutions of such problems were available. In 1980s- 
1990s, a considerable progress was achieved following the development 
of the theory of Whitham deformations [20] and the analysis of ma- 
trix Riemann-Hilbert problem representations of solutions of initial 
value problems, see [211 1221 123] and references therein. Most complete 
results,obtained by using this approach,were related to integrable equa- 
tions,for which linear operators from the associated Lax pair were self- 
adjoint and thus their spectrum was real. In [22].Bikbaev considered the 
case of the focusing nonlinear Schrodinger equation,which required the 
development of a much more complicated complex form of the theory 
of Whitham deformations. 

A completely rigorous approach for studying asymptotics of solutions 
of integrable nonlinear equations was introduced by Deift and Zhou 



4 J.XU, E.FAN, AND J.XU 

[9j(this approach was inspired by earlier works of Manakov [24j and 
Its [25];see |10] for a detailed historical review) and further extended 
by Deift,Venakides,and Zhou [2Sl U7\- This approach is based on the 
development of the nonlinear steepest descent method for Riemann- 
Hilbert problems associated with integrable nonlinear equations. Being 
originally introduced for studying initial value problems with decaying 
initial data, this approach was recently adapted by Buckingham and 
Venakides [2S] to problems with shock-type oscillating initial data for 
focusing nonlinear Schrodinger equation. A central role in this de- 
velopment is played by the so-called 5'— function mechanism allowing 
to deform the original Riemann-Hilbert problem to a form that can 
be asymptotically treated with the help of associated singular integral 
equations. 

The Riemann-Hilbert problem approach to initial value problems 
with nondecaying step-like initial data shares many issues with the 
adaptation of this approach for studying initial-boundary value prob- 
lems with non-decaying boundary data [221 EDI EI] -However, there is 
an important difference: in the latter case, the construction of the as- 
sociated Riemann-Hilbert problem normally requires the knowledge of 
spectral functions associated with overspecified initial and boundary 
data,which leads to the fact that results(in particular,the asymptotic 
results,see [29]) have, in a certain sense, a conditional character. As for 
the initial value problems of the type considered in this paper,the 
Riemann-Hilbert construction requires only initial data, and thus, the 
issue of overdetermination does not arise. 

In this paper,we consider a pure step-like initial value problem for 
the DNLS equation: 

iqt + qxx-iq^qx + -\q\^q = o, xeR,t>o, (1.6a) 

/ , , , f if X > 0, , ^ , 

q(x,0) = qo{x) = < ^.„ (1.6b) 

1 Ae-^'^^ if x<0, 
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where A > and i? G M are some constants. Kitaev and Vartanian 
got the leading order long-time asymptotic for the KN-type of DNLS 
equation with the decaying initial value,in |31], and the higher order 
long-time asymptotic in [36] . 

Since the DNLS equation (11 .Gap has a plane wave solution 

qP{x,t) = Ae-^'^^+^''^\ (1.7) 

with 

uj:=A^B-2B^ + —, (1.8) 

4 

which is consistent with fll.6b|) for x < 0,that is,g^(x,0) = qo{x),we 
assume that the solution q{x, t) of the initial value problem (ll.6ap eval- 
uated at any t > has the following behavior as x — ;■ ±oo: 

g(x,t) = o(l), x -)■ -l-cx), (1.9) 

q{x,t)=qP{x,t) + o{l), x^-oo, (1.10) 

where o(l) means sufficiently fast decay to O.This assumption can be 
justified a posteriori,by evaluating the large-x behavior of the solution 
of the Riemann-Hilbert problem formulated in Section 3. 

Recently, in [32],A.Boutet de Monvel,V.P.Kotlyarov, and D.Shepelsky 
considered the long-time dynamics of the initial value problem for the 
focusing nonlinear Schrodinger equation with step-like data. The strat- 
egy of the Riemann-Hilbert problem deformations that we adopt in 
this paper is similar,though not identical,to that in [28]. In particu- 
lar,the realization of the (7— function mechanism is different as well as 
the resulting asymptotic picture. 

As we have already mentioned, the main tool available now for study- 
ing rigorously the long-time asymptoitcs of solutions of initial and 
initial boundary value problems for integrable nonlinear equations is 
the asymptotic analysis of associated Riemann-Hilbert problems,whose 
construction involves dedicated solutions of the system of two linear 
equations,the Lax pair associated with the integrable nonlinear equa- 
tion. 
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For the DNLS equation f ll.6ap . a Lax pair is as follows p 

^^(x, t; k) = M{x, t; A;)^(x, t; k), 

(1.11) 
^i(x, t; k) = N{x, t] k)-^{x, t; k), 

where 

M(x, t; k) = -ik'^as + kQ + ^l^lVg, 

i 1 

N{x, t; k) = -2ik'^a3 + 2k^Q + 2A;^|g| V3 - ikQ^a^ + jlgl^cia + -{qqx - qqx)cr3, 

(1.12) 

with (Js = ( j , and \l/(x, t; fc) is a 2x2 matrix-value function,/!; G 

C is a spectral parameter, and the matrix coefficient Q is expressed in 
terms of a scalar function q: 



Q=\ : ], (1-13) 




It is well-known [3l] that this over-determined system of equations 
(II. lip is compatible if and only if q{x, t) solves the DNLS equation 
f lLHil) . 

In Section 2 we present these dedicated solutions (eigenfunctions) and 
associated spectral functions. All these functions are then used in Sec- 
tion 3 for constructing a basic Riemann-Hilbert problem,whose solu- 
tion gives the solution of the initial value problem (II. 6ap .f lL6b|) . Section 
4 develops the asymptotic analysis of this Riemann-Hilbert problem 
leading to asymptotic formulas for the solution of the original Cauchy 
problem (II. 6p . 

2. Eigenfunctions 

Let Q^ be defined by (I1.13P with q^ instead of q. A particular solution 
of the system (II. 11 p . with Q^ instead of Q,is given by 
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where 



with 



X{k) = d{k^-B-—y + k^A^, (2.2) 

n{k) = 2{k^ + B)X{k). (2.3) 

n^(k)-^ m + ^J 

The branch cut for X and </? is taken along the segment 

-fU^:={keC\kl-kl = B,kl<C^}, (2.6) 

where 7 = {A; G C|A;2-A;2 = B,kl < C^Jm^ > 0}, C^ = B + ^,ki = 
Kek and A;2 = IniA;. And the branches are fixed by the asymptotics: 

X(k) = P-B + 0(^), as A; ^00, 
k'^ 

Lp(k) = 1 + 0(— ), as A; — )■ 00. 
k 

We find that n{k) = 2k^ + u + 0(i),as A; ^ 00. We also find that 
ImX(A;) = is 

kik2{kl - kl - B) = 0, (2.7) 

which is on 

S:=MU2MU7U7. (2.8) 

Thus, for any t > 0, \&^(x, t; k) is bounded in x if and only if A; G S. 

Let q{x, t) be a solution of the Cauchy problem (ll.6ap .( lL6b| satisfy- 
ing the asymptotic conditions (ll.9p .( lLT0|) . and let Q{x,t) and Q^{x,t) 
be defined by (I1.13p . in terms of q and q^, respectively. Define the 2x2 
matrix-value functions ^j{x,t; k), j = 1,2, —00 <a;<oo,0<t<oo, 
as the solutions of the Volterra integral equations: 



Hi {x, t;k) = I+ I e'^ (^-^)"^ {kQfii){y, t; k)e-'^ (^-^)"=' , A;^ G M, 

) 

(2.9) 
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fX2{x,t;k) = e'("*-^")"^E(A;) (2.10) 

J ~oo 

where 

VP{x,y,t,k):=^^{x,t,k)[^P{y,t,k)]-\ 

Note that F^ can be written in the form 
where 



a + i{k'^-B-\)P -kA(5 

kA/3 a-i{e -B-^)(5 



G-{x^y^k)=\ --■\.: ^ .,,^, ^ ^: 



with 

r/ x,^/,M ^ sin[(-y — x)X(A;)l 

a = cos[(y - x)X{k)l (5 = ^J. ■ 

For any {x, y) G M.'^,G^{x, y, k) is an entire function of k with asymptotic 
behavior 

G'P(x,2/,fc) = e^(^-^)('='-^-^)"3[I + 0(^)], asA;^oo, ImA;^ = 0. 

rC 

The analytic properties of the 2x2 matrices fj^j{x, t;k), j = 1,2, that 
follow from (12 .Op and (12.101) are collected in the following proposi- 
tion. We denote by fij {x, t, k) and fij '{x, t, k) the columns oi fij{x, t; k). 

Proposition 2.1. The matrices fxi{x,t; k) and fX2{x,t] k) have the fol- 
lowing properties: 

(i) detfii{x,t,k) = fi2{x,t; k) = 1. 
(ii) The functions $(a;,t, k) and '${x,t, k) defined by 
^{x,t,k) := /ii(x,t,A;)e-^'^'^"3-2ifc4to3^ 
$(a;, t,k):= fi2{x, t; ;t)e-*-^W-3-itf^(fc)-3 _ 

satisfy the Lax pair equations U.ll\) . 
(iii) fj,i {x,t, k) is analytic in Im/c^ < and 

^i{\x,t,k)=[ \ + 0{-), ask ^oo, ImA;^ < 0. 
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(iv) fxl {x,t,k) is analytic in liak^ > and 

fif\x,t,k)=i ] +0{-), ask ^oo, ImA;^ > 0. 
\l J k 

(v) fj,2 (^) ^i ^) ^^ analytic in Imk^ > 0\'y, has a jump across 7, and 
fi'2\x,t,k) = I ] +0{-), as k^ 00, ImA;^ > 0. 

(vi) /i2 (x, t, k) is analytic in ImA;^ < 0\'y, has a jump across 7, and 
nf\x,t,k)={ ] +0{-), as k^ 00, Imfc^ < 0. 



(vii) Moreover, 



4'>(x.a^) = . + %^ + o(i, 



as k -^ 00 along curves transversal to the real and image axis, 
where 

(viii) t^2 {^, t, k) (k—E) 4 is boundary near k = E and ^2 (^) ^) ^) (^~ 
i?)* zs boundary near k = E. 

Since the eigenf unctions \l/(x,t, A;) and ^{x,t,k) satisfy both equa- 
tions of the Lax pair, we have 

^{x, t, k) = ^(x, t, k)S{k), k^ e M, (2.11) 

where S{k) is independent of {x,t). Since (see (I2.9p and (I2.10p for 
t = 0) 

^(x,0,A;) = e-*'='^'^^ for x > 0, 

$(x,0,A;) = e-*^""3^(fc)e-*"^('=)"^ for x < 0, 
we conclude that 

S{k) = ^'\0, 0, k)^{0, 0, fc) = $(0, 0, k) = E{k). (2.12) 
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Thus, we have 

S(k) = ( "^^^ ^^^^ ^ = ( "^^^ ^^^^ 1 (2 13) 

^ ' I -b{k) a{k) / I h{k) a{k) ' ' ^ ' ' 



where 



1 1 

a{k) = a{k) = -[(p{k) + 



b{k) = -m = iMk) ^ ^ 



(2.14) 



2'"^' ' ip{k)'' 
3. The basic Riemann-Hilbert problem 

The scattering relation (12. lip involving the eigenfunctions '${x,t, k) 
and $(x, t, k) can be rewritten in the form of conjugation of boundary 
values of a piecewise analytic matrix- value function on a contour in the 
complex fc— plane, namely: 

M+{x,t,k) = M_{x,t,k)J{x,t,k), A; G S, (3.1) 

where M±{x, t, k) denote the boundary vales of M(x, t, A;) according to 
a chosen orientation of S, and S = RUzRU7U7. 
Indeed, let us write (12.111) in the vector form: 

$W(x t, k) ^ ^(1) ^^ ^ r{k)¥^\x, t, k), 

^ ' (3.2) 



where 



^!-l^lhR = r{k)^^^\x,t,k) + ¥^\x,t,k), 
a[k) 



^W-^-A^I*^-«-f-^Wl- P'^) 



and define the matrix M(x, t, /c) as follows: 

( illl(|i^e^*^W ^(2)(^^^^^)g-ite(fc) )^ A; G {A; G C|ImA;2 > 0\7}, 

^^^'^'^^ " ^ ( ^«(x,t,A;)e^*^W ^^^^llf^e-**'^"'^ ), A; G {A; G C|ImA;2 < 0\7}, 

(3.4) 



where 



9{k) := 2k^ + -k\ (3.5) 
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Then the boundary values M+{x, t, k) and M_(a;, t, k) relative to S are 
related by ( 13.1 p . where 



J \X ^ Z^ rh j 



l-r\k) 
1 

1 



1 






(3.6) 



with 



(3.7) 



f{k):=r4k)-r_{k). 

The jump relation (13. ip considered together with the properties of 
the eigenf unctions listed in Proposition 1 suggests a way of representing 
the solution to the Cauchy problem (ll.6ap and ( ]1.6bp in terms of the 
solution of the Riemann-Hilbert problem, which is specified by the 
initial conditions (ll.6bp via the associated spectral function r{k). 

The solution q{x,t) of the initial value problem (ll.6ap and (ll.6bp 
can be expressed in terms of the solution of the basic Riemann-Hilbert 
problem as follows: 



q{x, t) = 2i lim {kM{x, t, k))i2- 



(3.^ 



where M is the solution of the following Riemann-Hilbert problem: 
Basic Riemann-Hilbert problem I. 

Given r{k),k'^ G R and f{k) = r+{k) - r_{k), A;^ G 7 U 7, and E = 
M U iR U 7 U 7, find a 2 x 2 matrix- value function M{x, t, k) such that 

(i) M{x,t, k) is analytic in A; G C\S. 
(ii) M{x, t, k) is bounded at the end points E and E. 
(iii) The boundary value M± {x, t, k) at S satisfy the jump condition 

M+{x,t,k) = M.{x,t,k)J{x,t,k), k e T, 

where the jump matrix J{x,t,k) is defined in terms of r{k) 
and f{k) by (M . 
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(iv) Behavior at oo 

M{x,t,k) =I + 0{-), ask-^oo. 

K 

If we try to analysis the long-time asymptotic behavior of the Gl-type 
of DNLS equation (11 .Gap and fll.6b|) with step-like initial value problem, 
this type of Riemann-Hilbert problem has a contradiction in the plane 
wave region. So we try to derive a new Riemann-Hilbert problem, which 
is similar to the type of nonlinear Schrodinger equation, to overcome 
this contradiction. That means we arrive at the following Riemann- 
Hilbert problem. 
We define 

N{x, t, k) = k-^M{x, t, k), (3.9) 

then the jump condition for A^ is 

N+{x,t,k) = N_{x,t,k)e-'^^'''+^^''^^UN{x,t,k). (3.10) 

introducing A = fc^ and control the branch of k as Signlmk = SignlmX, 
and define the modified scattering data p(A) = ^^, [T3] . 
Then 



X{X) = ^{X-B-^r + XA^ = y^(A -By + ^ + A^B, (3.11) 

n{X) = 2{X + B)X{X). (3.12) 

and the segment 

jUj:={XeC\Xi = B,Xl<D^}, (3.13) 

where j = {k e C|Ai = B,Xl < D^ImAs > 0}, D^ = A^B + ^, 
Ai = ReA and A2 = ImA. Let E = B + iD, then 7 = [E, B] and 
7 = [B, E]. And the jump condition for A^ is 

N+{x,t,X) = N_{x,t,X)e-''^^'^+^^'''^^UN{x,t,X). (3.14) 
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where 

Xp{X)e^m^) 1 

1 0\ 

Xf^X)e^mx) I j ' 

1 /(A)e-2**^W 
1 



JN{x,t, A) = < 



AG M, 

A G 7, (3.15) 

AG 7, 



where 



/(A) = p(A)+-p(A)_. 



(3.16) 



In other word,we have the following basic Riemann-Hilbert problem 

aA2 





d\--- 




•E 






1 


1 


• 






d 


7 


B 


Ai 


E = B + iD 
E^B-iD 


d\--- 




•E 





Figure 1 . The oriented contour E = M U 7 U 7. 

Basic Riemann-Hilbert problem II. 

Given p(A),A G M and /(A) = p(A)+ - p(A)_,A G 7 U 7, and S = 
M U 7 U 7, find a 2 X 2 matrix-value function N{x, t, A) such that 

(i) N{x, t, A) is analytic in A G C\S. 
(ii) N{x, t, A) is bounded at the end points E and E. 
(iii) The boundary value N±{x, t, A) at S satisfy the jump condition 

N+{x,t,X) = N_{x,t,X)JN{x,t,X), X G ^\{E,E,B}, 
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where the jump matrix J]^{x,t, k) is defined in terms of p(A) 
and /(A) by fl37[5|) . 
(iv) Behavior at oo 

N{x,t,X) =I + 0{-), as A ^ oo. 
A 



4. Long-time Asymptotics 

The representation of the solution q{x, t) of the initial value problem 
(11.61) in terms of the solution of an associated basic Riemann-Hilbert 
problem allows using the ideas of the asymptotic analysis of oscillating 
Riemann-Hilbert problems [HI EHl [IHl [HI 132] for studying the long-time 
asymptotics of q{x,t). The key fact leading to different asymptotics 
in different regions of the [x, t) half-plane is that the behavior of the 
jump matrix of the basic Riemann-Hilbert problem as a function of 
the large parameter t is different in these regions. Indeed, as seen on 
(I3.15p . this behavior is governed by the sign of Im^(A), which itself 
depends on (^ = ^. As we have already written, three regions are to be 
distinguished: 

(i) A Zakharov-Manakov region:^ > —B. 
(ii) A plane wave region:^ < —\/2D — B. 
(iii) An elliptic wave region:— v^D — B < ^ < —B. 

X = -At{B + ^/2D) ^ ^ x = -AtB,t 

■^ I 

^ plane wave > 
region ', 



elliptic wave ^x i 

region ^ n i 

\ 1 



Zakharov-Manakov 
region 



Figure 2. The different regions of the (x,f)— plane. 
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4.1. The Zakharov-Manakov region:^ > —B. In this region ^ > 
-B, we have Im^(A) > for aU A G 7 and Im6'(A) < for all A G 7. 
Therefore, the exponentials in the jump matrix J^, see fl3.15p . are 
decaying as t — )■ +00 for A G S\]R. 

This implies that one can follow the technique of asymptotic analysis 
proposed for the first time in [9]. The basic step of the procedure is 
a deformation of the original Riemann-Hilbert problem, with the help 
of the solution of an appropriate scalar Riemann-Hilbert problem, in 
order to obtain an equivalent Riemann-Hilbert problem whose jump 
matrix decays, in t, to a constant (in A) matrix. This leads to model 
Riemann-Hilbert problems whose solutions can be given explicitly. 

A particular feature of the Riemann-Hilbert problem under consid- 
eration is that the contour of the modified Riemann-Hilbert problem 
contains neither the real axis, where the jump matrix for the original 
Riemann-Hilbert problem oscillates with t, see fl3.15p . nor the finite 
parts 7 and 7. This happens due to the pure step-like initial condi- 
tions, which in turn implies that the associated spectral functions p(A) 
and Ap(A) can be analytically extended from the contour to the whole 
A-plane. 

4.1.1. First transformation. The first transform is as usual: 

M^\x,t,X) = N{x,t,X)5-'''{X), (4.1) 

where (|41]) 

is the solution of the following scalar Riemann-Hilbert problem: 

• 6{X) is analytic in C\(— oo,Ao], 

• 6{X) — 7- 1 as A — 7- 00, 

• 5(A) satisfies the jump relation 

5+iX) = 5_(A)(1 - Ap2(A)), A G (-00, Ao). (4.3) 
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Here, Aq is the stationary point of the phase function 6{X) = 2A^ + 4^A, 
that is, ^'(Ao) = 0: 

Then N^^^x,t, A) satisfies the jump condition 

N^^\x,t,X) = N^y{x,t,N)/ji\x,t,X), 
A G S(i) = S, 



(4.4) 



where 



that is 



tW 



Ji;>{x,t,X) = 61'JN5+''', 



tW 



Ji;>{x,t,X) = { 



^-ite&s I <5+ 



(l-Ap(A)2) -pM- 



_Xp_ 

5+5- 





5- 



5+5- ^ 5- 



'2itea:i 



5+ 



A gM, 

AG 7, 

A G 7. 
(4.5) 



From the Riemann-Hilbert problem of the (5, we can find 



r(i) 



-149(73 



-itOas 



1-Ap2 -p52 



1 



1 



J]v^(x,t,A) = < 



1-Ap2 53; 

1 

A/p2it6» 1 
1 /(52g-2^te 



^^^ 1-Ap2 







1 



A > Ao, 



A < Ao, 

AG 7, 
A G 7. 



(4.6) 



4.1.2. Second transformation. The next transformation is: 

N^^\x,t,X) = M^\x,t,X)G{X), (4.7) 
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where 



G(A) 



1 



l—Xp^ — 

1 



1 







Ap 1 „2Jte 1 
i-Ap2 53-e i 

1 -p(52g-2ite 

1 

1 



I, 



AgDi, 

XeD2, 

XeD,, 

AGD4, 

A G 1)5 U Z^e 



(4. 



The domains Di, . . . Dq are shown on the following Figure. 




Figure 3. The oriented contour S^^^ = Li U L2 U L3 U L4 



This new function A^*^^-* solves the equivalent Riemann-Hilbert prob- 
lem: 



N^+\x, t, A) = N^^\x, t, X)/n\x, t, A), 
AgS(2), 
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where 



4\x,t,X) = { 
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1 

1 

Ap 1 p2ite 1 

1 _p52g-2ite \ 

1 J 



1 







>^p„2ite -, 



AeLi, 
AeLa, 
AeLg, 
A G L4. 



(4.9) 



4.1.3. The last transformation. Now J]y (x, t, A) decays exponentially 
fast to the identity matrix, as t — )■ +00, and uniformly outside any 
neighborhood of A = Aq. Thus, we are in a situation where the asymp- 
totic analysis of ^41j works. Particularly, 

Ar(2)(a;^t,A) = Z(x, t, A)A^"'(x,t, A), 



where N°-'^{x,t,X) is a solution of the model problem explicitly given 
in terms of parabolic cylinder functions whereas Z{x, t, A) can be esti- 
mated: 

Z(x,t,A) = I + 0(^). 

Therefore, the final asymptotic result is as in [H] giving the main term 
of the asymptotic in terms of the modified reflection coefficient p(A): 

Theorem 4.1. (The Zakharov-Manakov region) In the region x > 
—AtB, the asymptotics, as t -^ +00, of the solution q{x,t) of the initial 
value problem 1^1. 6\) is described by the Zakharov-Manakov type formula 



g(x, t) = qas{x, t) + 0{——) 



(4.10) 
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where 



|a(Ao)P = ^ = -i^log(l-Ao|p(Ao)P), 

arga(Ao) = -3z/log2 - f + BigT{iv) - argr(Ao) + ^ /_'^ log |A - AoMlog(l - A|p(A)n, 

(4.11) 

4.2. The plane wave region: ^ < -\f2D - B. For x < -At{B + 
J2A^{B + ^)), that means, Im6'(A) is negative on 7 and positive on 7, 
which implies that the exponentials in ( I3.15P increase with t. Thus, the 
jump matrix Jjv for the Riemann-Hilbert problem does not converge 
to a reasonable limit as t — )■ 00. 

To bypass this difficulty, one deforms the Riemann-Hilbert problem 
in such a way that the phase Im6'(A) is replaced by another function, 
g{\)i providing suitable behavior of the modified jump matrix. The ex- 
tension of the nonlinear steepest descent method for Riemann-Hilbert 
problems, involving the (^-function mechanism was first proposed by 
Deift, Venakides, and Zhou, see [261 



4.2.1. The g function. A natural choice for a (^-function appropriate for 
the region adjacent to the half-axis x < 0, t = 0, is the phase appearing 
in the explicit expression for the eigenfunction \l/^, see (12. ip . associated 
with the potential q^. Setting 

g{x,t,\) = xX{\) + tVL{\), (4.12) 

where X{\) and VL{X) are defined in (13. lip and (I3.12p .we have 

The signature table for Im5f(A; ^) is the partition of the A-plane into 
maximal domains where the sign of Iing{X] ^) is constant. Its form can 
be controlled by the zeros of the differential dg{\). Indeed, 

MX) = 4 <^-''^^^^-^-> A (4.14) 
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where 






Thus, for ^ < -(5 + J2A^{B + ^)), ^± are both real. Moreover, 

B < ji_ < ji^ < —^. 

In what follows the signature table of the function hivg^X) for dif- 
ferent values of ^ plays a very important role. The lines of separation 
between the different domains are the real axile 

and the algebraic curve 

A^(Ai + = (Ai + 5 + 20[(Ai - B){\, + + ~^/^ ], (4.16) 
They are indeed given by Im5f(A) = 0. Because of 

Imr7(A)=4A2{(Ai + iJ + 20[(Ai-i?)(Ai+0 + ~^/^ ]-A^(Ai + 0} 
The equation fl4.16p can be written: 

A^(Ai + = (Ai + 5 + 20[(Ai - ^+)(Ai - M_)]. 

And the signature table of the function Im5'(A) is shown in the following 
Figure 4. 

The advantage of the signature table shown in Figure 4 is that 
there is a finite arc connecting the branch points E and E such that 
Imf7(A) = for all A along this arc. Since the jump matrix depends on 
t via exponentials of type e^'^^^^\ it is oscillatory along an arc where 
Im^(A) = 0. 

This suggests to deform the original contour 7 U 7 of the basic 
Riemann-Hilbert problem to a new contour 7^ U 7^ which depends on ^ 
and where lva.g{\) = 0, and to view ^(A), thus also ^^(A) as functions 
with branch cut 7^ U 7^. 
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Figure 4. The curves of Im^(A) = for x < -At{B + \ 2A^{B + ^)). 



Another important feature of g{X; ^) is that it has, up to a constant, 
the same large A asymptotic behavior as the phase function 0{X): 

g{X; = t{2\^ + 4eA + (7(00; 0) + O(^), A ^ 00, (4.17) 



where 



g{oo-i) = {u-ABi). 



(4.18) 



4.2.2. The first transformation. We put 

N^^\x,t,X) = e-^*^(~'«)"W(x,t,A)e-^(^^+2A2t-9(A)V3^ 

Then the matrix- value function A^*-^^ (x, t, A) satisfies the following Riemann- 
Hilbert problem: 



r(i) 



N^^\x,t,X) = N^_}\x,t,X)J)}>{x,t,X), AgSW =RU7<,U 



7<?, 
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with the jump matrix 



tW 



J'^>{x,t,X) = -{ 



1 - Xp\X) -p{X)e-^'aW 
Xp{X)e^'3W 1 

g-2J9-(A) Q 



A/(A) 







,2ig_(A) 

/(A) 

,2ig_(A) 



AG M, 
A £79, 



Ae%. 

(4.19) 
Here g±{X) are boundary values of (7 on 7^, U 7c,, and they are real. We 
also use the equation g+{X) = —g^{X). 

4.2.3. The second transformation. The next transformation is similar 
to the first transformation applied in the ZakharovManakov region, 
see Section 14.1.11 It involves the solution d{X) of the scalar Riemann- 
Hilbert problem 14. 3p but with //+ instead of Aq, where fi+ is the sta- 
tionary point of the new phase function g{X). With this new scalar 
function S{X), we set 

N^^\x,t,X) = A^W(x,t,A)r"^(A), 

Then the matrix- value function A^'-^-' (x, t, A) satisfies the following Riemann- 
Hilbert problem 

Nf\x,t,X) = N^y{x,t,X)/^\x,t,X), AeS(2) = S«, (4.20) 

(2) 

where J)^ (x^t.X) is defined as follows: 



jS\x,t,X) 



-igo-3 



1-Ap2 -pS^ 
1 



5^ 



-P_X2 



-«9o-3 



P 

1-Ap2"- 

T^'sl 1 " AP^ 

g2ig_(A) 



A > ^+, 

A < p+, 

Ae75, 
Ae%. 



(4.21) 
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4.2.4. The third transformation. The subsequent transformation 

N^^\x,t,\) = M^\x,t,X)G{X), 



involves G{X) defined similarly to (14. 8p . with t9 replaced by g and Aq 
replaced by /i+. Then N^^\x,t, A) satisfies the jump relation 



r(3)/ 



r(3)/ 



r(3)/ 



N^^'{x,t,X) = Nr{x,t,X)J}^'{x,t,X), 
across to the contour 



S(3) = Li U L2 U L3 U L4 U 7g U %, 



shown in Figure 5. 



i\ / 


+ 


/L4 




D, \E \/ 


/ 


^4 




si JM-Ai+\ / 


^ 








\ 






i?2 / \ 




\ i?3 




/^ \ 


- 


^\ 




L2 \ 


i^6 


iv3 





Figure 5. The contour S^^) = LiUL2UL3UL4U7gU7g of 
the Riemann-Hilbert problem for A^'-^^ for x < —At{B + 



A2^ 



2^2(5 + f)) 



And we notice that 

l.For A G L1UL2UL3UL4 the jump matrix j\^ [x, t, A) decays to the 
identity matrix, as t ^^ 00, exponentially fast and uniformly outside 
any neighborhood of A = /i+. 
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(3), 



2. For A G 7^, the jump matrix Jj^ {x,t, A) factorizes as 



1 (t^)-^'^-'^^-^'^ 

1 



e-2i9(A) Q 

A/(A)r2(A) e2^9-{A) 



1 (^)_,52,2.,_(A) 

1 

(4.22) 



3. For A G 7g, the jump matrix J)^ (x,t, A) factorizes as 
1 \ / e-^'9-W f{X)6\X) \ [ 



e2^^-W y 1^ (^)+5-2e2^5-W 1 



1 

(A) 

(4.23) 



4. Using the identities 



l + A/( 
l + /( 



-P 



1-Ap2' 
Ap , 



we find 



7(3) 



'J N l"^' '^1 ^) — * 



1-Ap2^" "' 

-{\f)-\\)6\\) 

\f{\)5-\\) 

/(A)<52(A) 

-ri(A)r2(A) 



-^^7^, 



^^79, 

(4.24) 
In order to arrive at a Riemann-Hilbert problem whose jump matrix 
does not depend on A, we introduce a factorization involving a scalar 
function F(A) to be defined; 



r(3) 



J'^\x,t,\) 



F-\\) 
F+(A) 



i \ / F_(A) 
i i V FZ^{\) 



(4.25) 
in such a way that the boundary values F±{X) of F{\) along the two 
sides of 7g U 7c, satisfy 



F_(A)F+(A) 



-i\f{\)5-\X) A G 79, 

zri(A)r2(A) A G %. 



(4.26) 
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Indeed, once f l4.25p is satisfied, one can absorb the diagonal factors into 
a new piecewise analytic function wliose jump across 7gU7g is only the 
constant middle factor in f l4.25p . 

Thus, we arrive at the following scalar Riemann-Hilbert problem: 
Scalar Riemann-Hilbert problem. 
Find a scalar function F{X) such that 

• -F(A) and F~^{X) are analytic in C\{7g U 7^}. 

• F{X) satisfies the jump relation: 

p ^.^ ^. ^ f -^A/(A)r2(A) = al\X)aZ\X)VX6-\X), X G 7,, 

"" [ if-\X)6-\X) = a^{X)a_{X)VX6-\X), X e %. 

(4.27) 

where the contour jg U jg is oriented from E to E, and 

• -F(A) is bounded at A = cxd. 

Introducing 

"(^) -{"-"""' Tcv^ (-«) 

then the jump relation fl4.27p transforms to 



r log^(A) . . logHjX) ^ i^^^i^^, AG7.U7,, 

^ X{X) J+ ^ X(A) ^ ] iHli^, AgR. 



log\/A(5-^(A) _ 

(4.29) 

L \ c im 

X(A) 



The Sokhotski-Plemelj formula shows that this last jump relation is 
satisfied by 

X(A) /■ \og^/s + log6-^{s,^) ds f log ab{s) ds 

(4.30) 
Then F(A) is defined in terms of H{X) by (I4.28p . At A = 00 we find 

where 

2^ JlaU% X+{S) Jr X{s) 
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with 

S{X, = exp — j^ y—^ dX , (4.32) 

Using the relation 1 — Ap^(A) = a^^(A), we find a simpler expression 
for 0(0: 

27r y^+ A (A) J^^U79 ^+(A) 

4.2.5. The fourth transformation. The factorization (14.251) suggests a 
fourth transformation 

N^'^\x, t, A) = F"^(oo, ON^'^\x, t, A)F-"«(A, 0, 

Then we have 

n'^\x, t, A) = N^^\x, t, X)/j^\x, t, A) 

For A G 7g U 7g the jump matrix Jj^ {x, t, A) is constant 

Jjy [X, t, A) = J]^ = 

\ 2 
l.For A G 7g U 7g the jump matrix J}^ {x, t, A) is constant: 

J^\x,t,X) = JT'=(^. ' 

\ I 

2. For A G L U L, the jump matrix Jpf {x, t, A) decays to the identity 

/^\x,t,X) = I + Oi^^). 

4.2.6. The final transformation. Finally, we can express A^*^^^ in the 
form 

iV(4) (^x, t, A) = m^'ix, t, X)N"'°'^{x, t, A) , 

where N'^°'^{x,t, A) solves the model problem: 

iV!"°^(x, t, A) = iVf "") (x, t, A) Jr", A G 7, U %, (4.33) 

with constant jump matrix 

\ i 
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and N^'^'^x, t, A) = I + 0(r 5). 

As for the model problem, since v^(A)_ = iLp{\)+ on 7^ U 7^, its 
solution can be given explicitly in terms of v^(A): 

4.2.7. Back to the original problem. Let N*{x,t, X), * = (l),(2),(3),(4),mod, 
denote the solution of the Riemann-Hilbert problem RH*, and let 

m*2(x,t)= lim(AM*(2:,t, A))i2, 

A— >oo 

Then, going back to the determination of q{x, t) in terms of the solution 
of the basic Riemann-Hilbert problem, we have 

q{x,t) =2im{x,t)i2 = 2ie'^''^^'^'^^m^^\x,t)i2 

=2ie2'3(~'«)m(3)(x,t)i2 + 0(r^) (4.34) 

=2ie''9^°^'^^m^^\x,t),2F'\oo,C) + Oit"^) 

Taking into account that g{oo, ^) = ut — 4Bx, 2ini"^°'^{x, t)i2 = A and 
F^^(oo,(^) = e~^*'^''^''we arrive at the following theorem: 

Theorem 4.2. (T'lane wave region^ In the region x < —At{B + 



2A^{B + —)),the asymptotics, as t ^ +cxd, of the solution q{x,t) of 
the initial value problem U.6]) takes the form of a plane wave: 

q{x, t) = ^e2*M-Sx-0(O) + 0(r 5), t -^ +CX). (4.35) 

Remark 4.3. // we let C, -^ +00, then /i+ — > +cxd, then (f){^) — > cf), 



with 0=2^/ u" x^(x) '^-^' ^^'^ '^^^'^ '^^^ above equation ( 4-61\ ) reduce to 



q{x,t) = Ae^"^^^^ ^^ '^\ this is correspondence to our initial condition 
up to a phase shift. 
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4.3. The elliptic region:-4t(i? + y/2D) < x < -AtB. For the limit 



case (^0 = ~{.B + \ 2A^{B + 4~))) we have /i+((^o) = Ai-('Co), see Figure 



7, whereas for ^ > —{B+\2A^{B + — )), yu^ and //_ become non-real, 
complex conjugated numbers. As a result, the (7- function mechanism 
with (/(A; ^) as in the plane wave region fails. This shows that there is 
a break in the qualitative picture of the asymptotic behavior at ^ = (^o- 



4.3.1. The new g- function. A suitable (^-function for ^ > —{B+\ 2A^{B -\- —)) 
can be obtained as follows. First, we need to introduce a new real sta- 
tionary point /i((^) which must be a zero of the new differential dg. On 
the other hand we have to preserve the asymptotic behavior of the g- 
function for large A. To do so we must change the denominator of the 
differential dg. Thus the new differential takes the form: 

_ „ (A-,K))(A-,_K))(A-,,K)) , 
^(X-E){\-E){\-d{0)(\-d{0) 

where /u(0)/^±(0! ^"^^ d{^),d{^) are to be determined. 

li H = d = d, then the new differential coincides with the previous 
one, that is dg = dg, which is expected to hold for the value ,^0 of ^ 
limiting the two adjacent asymptotic regions. 

Now we consider dg as an Abelian differential of the second kind 
with poles at oo± on the Riemann-Hilbert surface of 



uj{X) = y/(A-E)(A-^)(A-d(0)(A-J(0), 
with 

E = B + iD, d(0 = rfi (0 + ^^2(0 

The branch of the square root is fixed by the asymptotics on the upper 
sheet: 

uj{X) = X^ + 0{X), A ^00+. 

We choose on this Riemann surface a basis {a, b} of cycles as follows. 
The 6-cycle is a closed clock-wise oriented simple loop around the arc 
'jE,d joining E and d. The a-cycle starts on the upper sheet from the 



GI-DNLS LONG-TIME ASYMPTOTIC WITH STEP-LIKE INITIAL VALUE 29 

left side of the cut 7£,d, goes to the left side of the cut 7^3, proceeds 
to the lower sheet, and then returns to the starting point. 
We can also write the Abelian differential dg{X) in the form: 

dgiX) = ^ ^'^"'%"'^^"' dX, (4.37) 



and normalize it so that its a— period vanishes. This determines Cq: 

dX 

j{X) 



/i^(A3 + c,A^ + c,A)^ 



Co = 3 — G K. 

Jd uj{X) 

We also require that g{X) has the same large-A behavior as the original 
phase function 9{X): 

g{X) = 2Xh + AXx + 0{l), A ^ oo+. 

This condition implies 

C2 = C - B - di, 
Define ^(A) as the sum of two Abelian integrals: 

g(\i) = 2( r+ /j^l±i2^1±iilA±f£rfA. (4.38) 

Je Je i^l^J 

Then it evidently has real fe— period 

Now notice that ^(A) can be written as a single Abelian integral 

j,A).4r^l±£2^!±fA±f^dA 

Je ^l^J 

and indeed 

Bg= dg. 

The large-A asymptotics of g{X, ^) can now be specified as 



g{X,0 = ^XH + 4CXt + g{oo,0 + 0{X 



-i\ 
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where 

gioo, = t(2( r + r)[ ^' + ''^' + ''^ + '' -iX+0]dX+2D^-2B^-m) 

JE Je ^[^) 

(4.40) 
is a real function of ^. 

Remark 4.4. For ^ = —B, if we set fi{—B) = di{—B) = B and 
dii-B) = D, that is, d{-B) = E and d{-B) = E, then g{X,-B) 
coincide(up to a constant) with 6{X, —B): 

g{X,-B) = e{X,-B) + 2\E\'^. 

which provides matching at the interface with the Zakharov-Manakov 
region. 

In order to define /i, fi± and d as functions of ^, let us compare 
the forms fl4.36p and fl4.37p of the differential dg. This gives (/i± = 
/ii ± z/i2) : 

/i + 2/ii — di = B — ^, 

2^/xi + ^l + ^l + {^- B)d, - \dl = \D^ - Bi, 

1^{A + l4) = -coiC,di,d2). 

The local expansion of g{X) at A = c? is of the form 

^(A) = B^ + g^iX- df'^ + g2{X - df + ■■■, 

where Bg is real. The signature table for Img'(A) must have three 
branches of the curve Img'(A) = going out from the point d, see 
Figure 6. Indeed: 

• Since g{E) = 0, one branch should connect d with E. 

• There should exist a branch separating the basins of + and — 
near the real axis. 

• Since g{X) behaves like ^(A) for large A, there should be an in- 
finite branch going to infinity along the asymptotic line ReA = 
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Figure 6. The curves of Im5f(A) = for -4t(i? + 



2A^{B + ^)) <x< -AtB. 



Therefore, we arrive at the requirement gi = 0, that is 

(A - rf)^/V(A)|w = ,(^- MO)(^-/^-(0)(^-M O) ^ 

^ ' y^^iA-«f ^^X-E){X-E){d-d) 

The fact that /i is real imphes that fi+ = d and /i_ = d, which finally 
leads to the following ansatz for dg{X): 



where li{^), di{^) and ^2(0 {d = di + id2, ^2 > 0) satisfy the equations: 



H = B-^-di, (4.41a) 

dl = D'^-2{B-ij,){B-di), (4.41b) 

^ ' 2(A-/i)rfA = 0. 



,^_,^ y(A-5)2 + D2 
Recall that (14. 4 lap and (I4.41bl) follow from the requirement that 

dg{X) = {4X + A^ + 0{X-''^))dX, as A ^ oo. 



(4.41c) 



while (I4.41cp is the normalization condition J^ dg{X) = 0. 
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Substituting f l4.41ap and ( ]4.41bl) into f l4.41cp yields an equation re- 



lating implicitly di and ^. In terms of the variables u and v, where 

B-di ^ + B 



u 



D ' 2D 

this equation reads 






T{u,v)= I \l- '—^ -^ {ir + 2v-u)dr = 0. (4.42) 



which is considered for < v < ^ and m > 0. It is easy to check 
that J^iO^v) = 4ti(and thus J'{0,v) > for v > 0), J'{+oo,v) < 
0,J^(0,0) = J^(^,^) = and Mu,v) < for {u,v) ^ (^,^). 
Therefore, (I4.42p determines a unique function u = u{v),v G [0,^] 
such that u{0) = and u{^) = ^. Consequently, we have that 
the system fl4.4ip determines uniquely di{^),d2{^) and /i(0) such that 
di{-B - V2D) = B + ^/2D and di{-B) = B. 

We have now specified a (7— function ^^(A) whose signature table is as 
in Figure 8. Hence, we can begin deforming the basic Riemann-Hilbert 
problem. 

4.3.2. The first deformation. We deform the part 7 U 7 of the contour 
of the basic Riemann-Hilbert problem into a contour 'yE,E connecting 
E and E in such a way that it contains: 

(i) Two arcs 7^ and 7^ connecting, respectively, E with d and d 

and E, and where ling{X) = 0; 
(ii) An arc 7^ connecting d and d, passing through /i, and along 
which Img{X) < for ImA < and Im5f(A) > for ImA > 0. 

Supplying j^^^ = 7/^ U 7^ U 7^ with the orientation as going from E to 
E, we fix the branch of g{X) as having a jump across ■Je,e'- 



g{X)^ + g{X)_=0, 


A G 7d U 7d 


g{XU-g{X)^ = B^, 


A €7^, 


with ImB^ = 
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Figure 7. The contour S^^) = Li U L2 U L3 U L4 U 7^ U 



7d U 7^ U 7^ for -AtiB + J2A^{B + f )) < x < -AtB. 



4.3.3. The second transformation. The further series of transforma- 
tions 

N{x, t, A) ^ N^^^ (x, t, A) -^ N^^^ (x, t, A) ^ N^^^ {x, t, A) 

is similar to that for the plane wave region but 

(i) with g{X) replaced by g{\), 

(ii) with /i, which is the real stationary point of ^'(A) instead of /i+, 
(iii) with the partition into domains with boundaries L as shown 
in Figure 7. 

The jump matrix Jl^'{x,t, A) is as follows: 

• For A G Lj at a fixed positive distance from the stationary 
point A = /i((^), 



r(3) 



J}^'{x, t,X)=I + 0(6"^*) as t ^ +00 
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• For A G 7^ we have 



4\x,t,X) = { 



^-itB-g Q 



A/(A)r2(A)e^*(^+W+^-W) e^*^« / ' 



its, 



'« /(A)<52(A)e-^*(§+W+s-W) 







gitBg 



ImA > 0, 



ImA < 0, 



r(3) 



(4.43) 

Thus, away from d,fi and J and as t — )■ +oo, j]y (x,t, A) is 
close to a diagonal matrix: 

itBf, 





4'(a:,t,A) 



^itBg 



0{e 



-et\ 



t -^ +CX). (4.44) 



r(3)/ 



For A G 7dU7d, similarly to the plane wave region, J)^ {x,t, A) 
reduces to 



r(3) 



J'j^'{x,t,X) 



-f'^(X)5\X) 

A/(A)r2(A) 

f{X)S\X) 

-Xf-\X)6-\X) 



A G 7d, 



A G 7d, 

(4.45) 

In order to arrive at a Riemann-Hilbert problem with a jump matrix 
independent of A, we proceed as in the plane wave region. 

Scalar Riemann-Hilbert problem. We are looking for a scalar 
function -F(A) analytic in C\7d U 7^ such that 

F_(A)F+(A) = h{X)^/XS-\X), A G 7d U 7,, (4.46) 



where 



h{X) 



iVxfiX), 



Ae7s, 



(4.47) 



iVx' f-\k), Xe%. 

After solving this scalar problem, JV(x, t, A) can be factorized as in 
04.25p . This factorization allows absorbing the diagonal factors into a 
new Riemann-Hilbert problem with constant jump matrix on 7^ U 7^. 

However, an important difference with the plane wave region is that 
now the jump conditions (14.461) for -F(A) are specified on two disjoint 
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arcs. This implies that in order to arrive at a jump condition in additive 
form, we are led to use 

uiX) = ^(X-E){X-E){X-d{0){X-diO) 

Indeed, f l4.46p can be rewritten as 

JogF(A). rlogF(A). \ogh{X) 



^ u{X) '^ ' uj{X) ' w+(A) 

and thus for F{X), we have 



A G 7rf U 7rf, (4.48) 



F(A)=exp{^/ logM^ji.^ ^4_4g^ 



27ri ./^^U7d W-*-(s) s - X' 



But now F(A) has an essential singularity at infinity: 

F{X) = F^e'^\l + 0{X-^)), A^oo. 
where 



and 



with 

e. = '^+^ + " + '1, (4.51) 

To account for this singularity, let us introduce the normalized, that 
is, its a— period vanishes, Abelian integral it;(A) of the second kind with 
simple poles at oo-t: 

w[X) = / r az, 



where ei is the same as in fl4.5ip and cq is determined by the condition 
JJwiX) = 0: 



eo 



(■d dz 

■id OJrz) 

The large-A expansion of w{X) is of the form 

w{X) = X + Woo{0 + O{X'^), A^cx), 
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where 

Woo = / ^ 1 02; - E 

^' "'-'' , (4.52) 

The jump conditions for w{X) are as follows: 

w+{X) + w^{\) = 0, AG7dU7rf, 
w+{X) - w^{\) = By,, A G 7;,. 
Here B^, is the &— period of to (A): 

Jh Je ^(Z) Je Je ^(Z) 

(4.53) 



6 Je ^{Z) Je Je ^{Z} 

Now introduce 



F(A) = F(A)e-*^"'(^), (4.54) 

This new function is clearly bounded at A = 00: 

F(cx),0 = e^^(«\ (4.55) 

with 

0(0 = ^ I {s-e,) log [h{s)5-\s, 0]^ - A(0^oo(0- 

Also, F{\) has the same jumps as F{\) across 7^ and 7^. On the other 
hand, the price for introducing the exponential factor in 04.54p is that 
F(A) has a jump across 7^: 

F_(A) ' ^'^^ 

Now we can absorb F{\) into the Riemann-Hilbert problem for A^^^-* (x, t, A): 

A(^) (x, t, A) = F'^^ (cx)) A^^^^) (x, t, A)F-"« ( A) , 

which leads to the jump conditions 

N^^\x, t, A) = N^_^\x, t, X)/^\x, t, A), 
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where 



J^^\x,t,\) 






Xe LUL. 



with 



J. 



{mod) 



N 



1 




1 oj- 




^-itBg-iABn, 








^UBg+iAB 



A G 7d U 7d, 



AG 7^, 



(4.56) 



4.3.4. The model problem. Thus, we arrive at the model Riemann- 
Hilbert problem: 



jmodi 



jmod 



iV^°^(a;,t, A) = iV!"°'^(a;,t, A)Jr'(a;,t, A), A G 7dU7dU7^, (4.57a) 



N"'°''{x,t,X)=I + 0{-), A^oo. 

A 



(4.57b) 



The solution of this model Riemann-Hilbert problem approximates 

A^(4)(x,t,A): 



N^^\x,t,X) = {I + 0{t-'^))N"'"\x,t,X), 



(4.58) 



The model problem (14.571) can be solved in terms of elliptic theta func- 
tions. Let 



[/(A) = - 



dz 



cJe ^{z) 
be the normalized Abelian integral, that is 

'■^ dz 



c = 2 



d ^W 



Then, define 



r = r(0 = - 



dz 



c Je uj[z) 
with Imr > 0. Furthermore, the following relations are valid: 

f/+(A)-Ff/_(A)=0, AG7d, 

f/+(A)-ft/_(A) = -l, AG7d, 
f/+(A)-f/_(A) = r, AG 7m, 



(4.59) 



(4.60) 
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Next, define 

where the branch is fixed by specifying the branch cut 7^ ^ ^-nd the 
behavior as A — )■ 00; 

zy(A) = 1 + ^^^ + 0(A-2), A ^ 00. 

Along the cut, we have 

1 -'^-(A), AG 7^. 
Finally, introduce the theta function 

mail 

and define the 2x2 matrix-value function 0(A) = 0(t,(^,A) with en- 
tries: 

1 ^6,[U(X)-Uo~\-^-^] 



■Kirm^ +2nimz 



01l(A) - |[//(A) + ^j 03[t/(A)-C7o] 

<Ji2(Aj - ^[p[A) - ^ e^iuw+Uo 



P) (\)-^\iA\) 1 ^ 03[UW+Uo-^-^-^] 

<J2UAj - ^[l^[A) - ^jp^J g^iu{X)+Uo] ' 

<J22lAj - 2i/^lAj + ^;p^J g.^lu(X)-Uo] ' 

where Uq is to be chosen so that the unique zero of 93{U{X) — Uq), as a 
function on the Riemann surface, lying on the first sheet is compensated 
by the zero of z/(A) + j]^ where O^^U^X) + ^0) has no zero on this sheet. 
Setting 

Uo = UiEo) + 1 + ^, 

where 

Ed -Ed 

En 



E-E+d-d 
satisfies this requirement, and thus 0(A) can be viewed as a function 

analytic in C\7^ £;. On the other hand, due to the properties of theta 

function: 

e,{-z) = e^{z), e^{z + 1) = d,{z), d,{z ± r) = e--^^^^-^^e,{z) 
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0(A) satisfies the jump conditions f l4.57ap -( l4.56p of the model Riemann- 
Hilbert problem. Taking into account the normalization condition 
(I4.57bp . the solution of the model Riemann-Hilbert problem is given 
by 

4.3.5. Back to the original problem. Now, following the sequence of 
equations of type fl4.34p (with g and F replaced, respectively, by g and 
F) and taking into account the equations g and F, and the explicit 
formula for n'^'^{x, t, X) 

and -F~^(oo) = e~'^^'^^^\ we obtain the asymptotics in the region —At(B+ 
2A^{B+^)) <x< -AtB. 



A^' 



Theorem 4.5. (T^lliptic wave region^ In the region —At{B+J2A^{B + ^)) < 
X < —AtB, the asymptotics, as t ^ +oo, of the solution q{x,t) of the 
initial value problem U.6]) takes the form of a modulated elliptic wave: 



M¥ + ^ + ^'^(O]0s[v4O '' 



(4.61) 
Here Bg,Bw and A are functions of the variable ^ = fj: defined, 
respectively, by i4.39\j , ^.53\ ) and ( [^.50| j, and V±{^) = f/o + | ±f/(oo). 
Furthermore, 

zez 
is the theta function of invariant r = r(^) defined in ^-59 ), 



oo. 



and the phase shift 0(^) is given by 

, ... ^J_ f [s - ei(0 - uj^m hg[his)^sS-\s,0] 
^ ' 27r 4,^^_ [(. - E){s - E){s - d{0){s - mW^^ 
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where 



a+(A)a_(A), A G 7j 
5(A,0=exp{^/_-«)Mi±Amd.}. 



and ei{^),uJoo and fi{^) are defined, respectively, by ( [-^.5i[ j, U-^^ and 
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